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This article derives an optimal (i.e., unbiased, minimum variance) estimator for the pseudo- 
detector strain for a pair of co-located gravitational wave interferometers (such as the pair of LIGO 
interferometers at its Hanford Observatory) , allowing for possible instrumental correlations between 
the two detectors. The technique is robust and does not involve any assumptions or approximations 
regarding the relative strength of gravitational wave signals in the Hanford pair with respect to 
other sources of correlated instrumental or environmental noise. 

An expression is given for the effective power spectral density of the combined noise in the pseudo- 
detector. This can then be introduced into the standard optimal Wiener filter used to cross-correlate 
detector data streams in order to obtain an optimal estimate of the stochastic gravitational wave 
background. In addition, a dual to the optimal estimate of strain is derived. This dual is constructed 
to contain no gravitational wave signature and can thus be used as an "off-source" measurement to 
test algorithms used in the "on-source" observation. 

PACS numbers: 04.80.Nn, 04.30.Db, 95.55.Ym, 07.05.Kf, 02.50.Ey, 02.50.Fz, 98.70.Vc 



I. INTRODUCTION 

Over the past few years a number of long-basclinc 
interferometric gravitational wave detectors have begun 
operation. These include the Laser Interferometer Grav- 
itational Wave Observatory (LIGO)detectors located in 
Hanford, WA and Livingston, LA Q; the GEO-600 de- 
tector near Hannover, Germany the VIRGO detector 
near Pisa, Italy Ja|; and the Japanese TAMA-300 de- 
tector in Tokyo j4j. For the foreseeable future all these 
instruments will be looking for gravitational wave signals 
that are expected to be at the very limits of their sensi- 
tivities. All the collaborations have been developing data 
analysis techniques designed to extract weak signals from 
the detector noise. Coincidences among multiple detec- 
tors will be critical in establishing the first detections. 

In particular, LIGO Laboratory operates two co- 
located detectors sharing a common vacuum envelope at 
its Hanford, WA Observatory (LHO). One of the two de- 
tectors has 4 km long arms and is denoted HI; the other, 
with 2 km long arms, is denoted H2. This pair is unique 
among all the other kilometer-scale interferometers in 
the world because their co-location guarantees simulta- 
neous and essentially identical responses to gravitational 
waves. This fact can provide a powerful discrimination 
tool for sifting true signals from detector noise. At the 
same time, however, the co-location of the detectors can 
allow for a greater level of correlated instrumental noise, 
complicating the analysis for gravitational waves. 

Indeed, it may not be feasible to ever detect a stochas- 
tic gravitational wave background, or even establish a 



significant upper limit, via cross-correlation of HI and 
H2, due to the potential of instrumental correlations. 
However, even though it may not be profitable to cor- 
relate these co-located detectors, the data from HI and 
H2 should be optimally combined for a correlation anal- 
ysis with a geographically separated third detector (such 
as LI, the LIGO Livingston detector). 

For the H1-H2 detector pair, properly combining the 
two data streams will always result in a pseudo-strain 
channel that is quieter than the less noisy detector. In 
the limit of completely correlated noise, this combina- 
tion could, in principal, lead to a noiseless estimate of 
gravitational wave strain. In the other limit where the 
detector noise is completely uncorrelatcd, the two detec- 
tor outputs can of course be treated independently and 
combined at the end of the analysis to produce a more 
precise measurement than either separately, as done in 
Section V.C. of Ref. It is the more general inter- 

mediate case, where there is partial correlation of the 
detector noise, that is the subject of this paper. 

We show that it is possible to derive an optimal — 
i.e., unbiased, minimum variance — strain estimator by 
combining the two co-located interferometer outputs into 
a single, pseudo-detector estimate of the gravitational 
wave strain from the observatory. An expression is given 
for the effective power spectral density of the combined 
noise in the pseudo-detector. This is then introduced into 
the standard optimal Wiener filter used to cross-correlate 
detector data streams in order to obtain an estimate of 
the stochastic gravitational wave background. 

Once the optimal estimator is found, one can subtract 
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this quantity from the individual interferometer strain 
channels, producing a pair of null residual channels for 
the gravitational wave signature. The covariance matrix 
for these two null channels is Hermitian; it thus pos- 
sesses two real eigenvalues and can be diagonalized by 
a unitary transformation (rotation). Because the covari- 
ance matrix is generated from a single vector, only one 
of the eigenvalues is non-zero. The corresponding eigen- 
vector gives a single null channel that can be used as an 
"off-source" channel, which can be processed in the same 
manner as the optimal estimator of gravitational wave 
strain. 

The technique described here is possible for the pair of 
Hanford detectors because, to high accuracy, the grav- 
itational wave signature is guaranteed to be identical 
in both instruments, and because we can identify spe- 
cific correlations as being of instrumental origin. Coher- 
ent, time-domain mixing of the two interferometer strain 
channels can thus be used to optimal advantage to pro- 
vide the best possible estimate of the gravitational wave 
strain, and to provide a null channel with which any grav- 
itational wave analysis can be calibrated for backgrounds. 

The focus of this paper is the development of this 
technique and its application to the search for stochastic 
gravitational waves. However, it appears that any other 
search can exploit this approach. 

In Section[n]we discuss the experimental findings dur- 
ing recent LIGO science runs which motivated this work 
to extend the optimal filter formalism in the case where 
instrumental or environmental backgrounds are corre- 
lated among detectors. In Section ITTT1 we introduce the 
optimal estimate of strain for the pair of co-located Han- 
ford interferometers. In Section IIVI we then introduce 
the dual null channel. Then in Section[V]we apply these 
formalisms to measurement of a stochastic background 
and consider limiting cases that provide insight to under- 
standing the concept. Finally in Section IvH we discuss 
the implications of these results and estimate the effects 
of imperfect knowledge of calibrations on the technique. 
Appendices ^ contain derivations of formulae used in 
Sec. E| 



II. INSTRUMENTAL CORRELATIONS 

Early operation at LIGO's Hanford observatory has 
revealed that the two LHO detectors can exhibit instru- 
mental cross-correlations of both narrowband and broad- 
band nature. Narrowband correlations are found, e.g., at 
the 60 Hz mains line frequency and harmonics, and at fre- 
quencies corresponding to clocks or timing signals com- 
mon in the two detectors; these discrete frequencies can 
be identified and removed from the broadband analysis of 
a stochastic background search, as described in Ref. 0. 
Broadband instrumental correlations, on the other hand, 
are more pernicious to a stochastic background analysis; 
the following types of relatively broadband correlations 
have been seen at LHO: 



• Low-frequency seismic excitation of the interfer- 
ometer components, up to approximately 15 Hz; 
at higher frequencies, the seismic vibrations are 
not only greatly attenuated by the detectors' iso- 
lation systems, but they also become uncorrelated 
over the distances separating the two interferome- 
ters. These correlations are not directly problem- 
atic, since they arc below the detection band's lower 
frequency of 40 Hz. 

• Acoustic vibrations of the output beam detection 
systems. 

• Upconvcrsion of seismic noise into the detector 
band: intcrmodulation between the mains line fre- 
quencies and the low-frequency seismic noise pro- 
duces sidebands around the {60 Hz, 120 Hz, . . . } 
lines that are correlated between the two detectors. 

Magnetic field coupling to the detectors is another po- 
tential source of correlated noise, though this has not yet 
been seen to be significant. 

The analysis of the first LIGO science data (SI) for a 
stochastic gravitational wave background showed sub- 
stantial cross-correlated noise between the two Hanford 
interferometers (HI and H2), due to the above sources. 
This observation led to disregarding the H1-H2 cross- 
correlation measurement as an estimate of the stochastic 
background signal strength. Two separate upper limits 
were obtained for the two transcontinental pairs, Ll-Hl 
and L1-H2 (LI denotes the 4 km LIGO interferometer in 
Livingston, LA). These were not combined because of the 
known cross-correlation contaminating the H1-H2 pair. 

Here, we show how to take into account such local 
instrumental correlations in an optimal fashion by first 
combining the two local interferometer strain channels 
into a single, pseudo-detector estimate of the gravita- 
tional wave strain from the Hanford site, and then cross- 
correlating this pseudo-detector channel with the single 
Livingston detector output. In doing this, we obtain a 
self-consistent utilisation of the three measurements to 
obtain a single estimate of the stochastic background 
signal strength fi gw . In order for this to be valid, the 
reasonable assumption is made that there are no broad- 
band transcontinental (i.e., Ll-Hl, L1-H2) instrumental 
or environmental correlations. This has been empirically 
observed to be the case for the SI, S2 and S3 science 
runs when the Ll-Hl and L1-H2 coherences are calcu- 
lated over long periods of time (the SI findings are dis- 
cussed in Q; S2 and S3 analyses are still in progress at 
the time of this writing). 

It is important to point out that the technique pre- 
sented here is robust and does not involve any assump- 
tions or approximations regarding the relative strength 
of gravitational wave signals in the H1-H2 pair with re- 
spect to other sources of correlated instrumental or envi- 
ronmental noise. Since SI, the sources of environmental 
correlation between the Hanford pair have been largely 
reduced or eliminated. However, as the overall detector 
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noise is also reduced, smaller cross-correlations become 
significant, so it remains important to be able to opti- 
mally exploit the potential sensitivity provided by this 
unique pair of co-located detectors. 



III. OPTIMAL ESTIMATE OF STRAIN FOR 
THE TWO HANFORD DETECTORS 

Assume that the detectors HI and H2 produce data 
streams 

s Hl {t) := h(t)+n Hl (t), (3.1) 

s H2 (t) := h(t)+n H2 (t), (3.2) 

respectively, where h(t) is the gravitational wave strain 
common to both the detectors. In the Fourier domain, 

s Hl (f) = Hf)+n Hl (f), (3.3) 
SH 2 (f) = h(f) + n Hz (f), (3.4) 

where we defined the Fourier transform of a time do- 
main function, a(t), as a(f) := f_ co dte~ l2 ' n 'f t a(t). Also 
assume that the processes generating h, tlh 11 nn 2 are 
stochastic with the following statistical properties: 



nnAf)) = (h(f)) 


= o, 


(3.5) 




= o, 


(3.6) 


(n* Hi {f)n Hi {f)) 


= PZ iH ,(f)Hf-f), 


(3.7) 


(h*(f)h(f)) 


= Pa(f)S(f-f), 


(3.8) 


<3*r 4 (/)**,(/')} 


■■= P HiHj {f)S{f-f) 


(3.9) 




= (P£ iHj (f) + Pn(f)) 


(3.10) 




xS(f-f) 


(3.11) 




■= PhAD, 


(3.12) 


Ph.hAI) 


■■= P Hi {f), 


(3.13) 




P Ht H 3 (f) 

VPhADPhAJ) ' 


(3.14) 




:= |PW/)| 3 , 


(3.15) 


Pn(f) 


« PH t (f), 


(3.16) 



where i = 1,2 and the angular brackets (...) denote en- 
semble or statistical averages of random processes. Note 
that Eqs. (|3.9|) and 13.13f) signify the measurable cross- 
power and power spectra while Eqs. (|3.7() and l|3.12|l refer 
to intrinsic noise quantities that cannot, in principle, be 
isolated in a measurement. Often, Eq. (|3.16|) is assumed 
in order to identify instrument noise power with the mea- 
sured quantity. Note also that the coherence PHiHj(f) is 
a complex quantity of magnitude less than or equal to 
unity, and that P H] hA!) = P^H^f)- 

Now construct an unbiased linear combination of 

Mf) :=a(f)s Hl (f) + (l-a(f))s H2 (f). (3.17) 



If Sh{/) is also to be a minimum variance estimator, 
where 



Var( Sff ) := (T H (f)s H (f)) = P H (f) S(f - /') , (3.18) 
with 

PhU) - |5(/)| 2 P5 1 (/) + |l-5(/)| 2 P5,(/) + 
+a*(f)(l-a(f))P^ lH Jf) + 

+5(/)(l-S*(/))P£ Ha (/)+ifc(/) > 

(3.19) 

then a(f) must have the following form: 

~ m = P H2 {f) - P Hl H 2 {f) 

P Hl (/) + P H2 (/) - (P Hl H 2 (/) + P* Hl H 2 (/)) ^ 

(3.20) 

The corresponding power of the pseudo-detector signal is 



H(J> P Hl (/) + P H2 (f) - (P Hl H 2 (f) + P* Hl H 2 U)) ' 

(3.21) 

It is important to note that the above expressions for 
a(f) and Pu{f) do not require any assumption on the 
relative strength of the cross-correlated stochastic sig- 
nal to the instrumental or environmental cross-correlated 
noise. In particular, the stochastic signal power Pq enters 
Phii Ph 2 i and PhiH 2 in exactly the same way, canceling 
out in Eq. (|3.20|l . implying that the above solution for 
a is independent of the relative strength of the stochas- 
tic signal to other sources of cross-correlated noise. In 
addition, Eqs. H3.20|l . (j3.21(l involve only experimentally 
measurable power spectra and cross-spectra (and not the 
intrinsic noise spectra), indicating that this procedure 
can be carried out in practice. 

Figure ^ shows plots of the strain spectral densities 
for Sh(/), and Sh 2 (/), representative of the SI 

data. The strain spectral density \Sh{J)\ is calculated 
from Eqs. (|3~T7|) and (|3~3U|) for both T Hl H 2 (f) = (i.e., 
an artificial case that assumes no coherence), and for the 
coherence ^H t H 2 {f) that was actually measured over the 
whole SI data run (see Fig- ■ The plots in Fig. ^sug- 
gest that the observed level of coherence during the SI 
run, F ~ 10 -5 , might be sufficiently low that one can 
simply combine the Ll-Hl, L1-H2 cross-correlation mea- 
surements under the assumption of zero cross-correlated 
noise (c.f. Eq. I|5.20|l ). The formalism developed in this 
paper allows a quantitative assessment of the effect of in- 
strumental or environmental correlations on combining 
independently analyzed results ex post facto. 
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Representative El strain sensitivities for 
s H = a ( f ) s H1 (f) + (1 - a ( f ) ) s H2 ( f ) 
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FIG. 1: Strain spectral densities (i.e., absolute value) of sii(f) (gray or dotted), SH 1 (f) (black), and 5h 2 (/) (dashed), represen- 
tative of the SI data. Top Panel: overlay of the individual spectral densities with that of the strain spectral density |Sh(/)| 
calculated with the SI run-averaged coherence, ^H 1 H 2 (f), an d with ^H 1 H 2 (f) — 0. On this scale, the left hand panel shows no 
discernible difference between the spectra for ^H 1 B 2 (f)> an d with ^H 1 H 2 (f) = 0, suggesting that even the level of coherence 
seen during the SI run might be sufficiently low to allow one to simply combine the Ll-Hl and L2-H2 cross-correlation mea- 
surements under the assumption of zero cross-correlated noise. The optimality of the estimate sjj (/) is visible here because it is 
always less than the smaller of Sfr 1 (/) or Sif 2 (/). The inset shows a blow-up of the region near one of the spectral features. On 
this scale the individual spectra can be discerned. Bottom panel: plot of the ratio of amplitude spectra for |Sh(/)| calculated 
with r_r/ 1 H 2 (/) as measured during SI and ^H 1 H 2 (,f) = (i.e., assuming no coherence). The difference between the two is very 
small except for the very lowest frequencies and at narrow line features. 



A. Limiting cases II. If Pjj 1 {f) = Pn 2 (f), then 

I. If PiW/)=0, then 

a( f) PhAI) (3 22) 5 (/) = o { ^ m\ ■ < 3 - 25 ) 

a ^> - P Hl (f) + Ph 2 U) ' [ ' Z-{p HlH2 (f)+p HiH2 (f)) 

SHif) = P Hl (f) + PH 2 (f) ' (3 ' 23) 

PHif) = pWBSty (3 - 24) 
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Run averaged H1-H2 coherence, T(f) , for SI 
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FIG. 2: H1-H2 coherence averaged over the whole SI data run. Note the substantial broadband coherence below 250 Hz and 
between 400 and 475 Hz. Low frequency seismic noise and acoustic coupling between the input electro-optics systems are 
considered to be the prime sources of this cross-correlated noise 0, • 



III. II p„,„,(f ) = 1 and J>„,(/) = P H ,(/), then 

W) = n 1 )',".. F "f r^m - p "- (f) ■ <3 - 26) 

IV. If P H2 (/) = 4P ffl (/) (which is the limiting design 
performance for HI and H2 due to the 2 : 1 arm length 
ratio), then 



2(2 -p HlH2 (f)) 



5-2( PHlH2 (f) + p* HiH2 (f)) 



(3.27) 



Note for this case that if the noise were either completely 
correlated (pH 1 H 2 (f) = 1 = 2) or anti-correlated 

(pH 1 H 2 {.f) = — 1 => = 2/3), then one could exactly 
cancel the noise from the signals . If the noise is 
uncorrelated (pH 1 H 2 (f) = => a(f) = 4/5), then the 
weighting of the signals from the two interferometers is 
in the ratio 4 : 1, as expected. 



IV. A DUAL TO THE OPTIMAL ESTIMATE OF 
STRAIN THAT CANCELS THE 
GRAVITATIONAL WAVE SIGNATURE 

In the previous section, an optimal estimator of the 
gravitational wave strain h was derived by appropriately 
combining the outputs of the two Hanford detectors. It 
is also possible to form a dual to this optimal estimate 



(denoted that explicitly cancels the gravitational 

wave signature. 

Starting with Eqs. (|3.3|) . (|3.4|) . and the optimal esti- 
mate ?£?(/), we construct the /i-subtracted residuals: 

Z Hl (f) ■= 8Hdf)-*H{f), (4- 1 ) 

ZH 2 {f) := s H2 (f)-s H (f). (4.2) 

Both z Hl (/), z H2 (f) arc proportional to n Hl (f)—nn a (/), 
although with different frequency-dependent weighting 
functions: 

z Hl [f) = (1 - S(/)) (n Hl (/) - n H2 (/)) , (4.3) 
zn 2 (f) = -a(f) (n Hl (f)-n H2 (f)) ■ (4.4) 

Figure [21 shows schematically the geometrical relation- 
ships of the signal vectors s# 4 (/) and 2h 4 (/). Once the 
best estimate sji(f) is subtracted from the signals, the 
residuals lie in the nH 1 -fiH 2 plane. (Here fiu 1 and fiH 2 are 
unit vectors pointing in directions corresponding to un- 
correlated detector noise.) Their covariancc matrix can 
then be diagonalized without affecting the gravitational 
wave signature h contained in s//(/). 
Now consider the covariance matrix 

HCally 5(f- /') := (z* Hi (f)z Hj (f)) ■ (4.5) 
Then one can show that 



Schematic of signals showing 
noise and pW components 
h 
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FIG. 3: Schematic showing how the HI and H2 signals may be represented in a 3-dimensional space of noise components for 
the two detectors and the common gravitational wave strain: {nn 1 , uh 2 , h}. The signals Ish^ (/) and 1h 2 (/) are not, in general, 
orthogonal if the coherence between the noise, nn 1 (f) and UH 2 {f), is non-zero. sh(/) is the minimum variance estimate of 
h(f) derived from Sif 1 (/) and sh 2 (/). Using 3h(/) as the best estimate of h(f), this signal can be subtracted from Sfr 1 (/) and 
SH 2 {f) t° produce the vectors ZHxif), ZH 2 (f) that lie in the fiH 1 -nH 2 plane. These vectors give rise to the covariance matrix 
||C Z (/)||. ZHiif) and zh 2 {/) are colinear and thus one of the eigenvectors of ||C Z (/)|| will be zero. The other corresponds 
to the dual of Sh{/), denoted zn{f), which is orthogonal to 3h(/), as shown in the figure. Note that it is necessary to first 
subtract the contribution of h(f) from the signals before forming the covariance matrix. 



|C,(/|| *(/-/') = 



(^ 1 (/)^ 1 (/')> (z* Hl (f)^H 2 (f)) 
(z*H 2 (f)z Hl (n) (zH 2 (f)zH 2 (f')) 

\l-a(f)\ 2 -5(/) + |5(/)| 2 
-5*(/) + |5(/)| 2 |5(/)| 2 

|l-5(/)| 2 -5(/) + |5(/)| 2 
-5*(/) + |5(/)| 2 |5(/)| 2 

|l-5(/)| 2 -5(/) + |S(/)| 2 



(4.6) 

((n* Hl (/) - n* H2 (f))(n Hl (/') - n H2 (/'))> (4.7) 
((^(/) -^ 2 (/))(^(/') ~sh 2 (/'))) (4.8) 



-S*(/) + |3(/)f 



X (P Hl (f) + PH 2 (f) - (P Hl H 2 {f) + P* Hl H 2 (f))) S(f - /') . 

Diagonalization of ||C Z (/)|| gives the eigenvalues: 
Ai = 0, 

A 2 = (P Hl (f) + PhM) - (P Hl H 2 (f) + P* HlH2 (f))) (l-5(/)-5*(/) + 2|5(/)| 2 ) 
The non-trivial solution corresponds to the desired "zero" pseudo-detector channel: 

v / l-5(/)-5*(/) + 2|5(/)| 2 , 



(4.9) 



(4.10) 
(4.11) 



(4.12) 
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where a{f) is given as before (c.f. Eq. I |3.20| l1. The power spectrum P z (f) of zii(f) is given by the eigenvalue A2 
above. , 



Figure 0| shows plots of the strain spectral densities for 
SfTi(/), and Sjj 2 (/), representative of the SI data, 
similar to Fig. ^ 



A. Limiting case for zero cross-correlated noise 

In the limit that the two detectors are uncorrelated 
(i.e., PH 1 H 2 (f) = 0), the expression for a(f) simplifies 
considerably (c.f. Eq. I|3.22)l ). In this limit, zh(J) and 
P z {f) become 



p,(.f) = 



(s H Af)-SH 2 (f)) 



P Hl (f) + PhAI) ' 

(4.13) 
(4.14) 



PhAD + PhM) ■ 

In particular, note that P z {f) satisfies the inequality 

max{P Hl (/), P H2 (/)} - min{P Hl (/), P H2 (/)} 

< Pz(f) < max{P Hl (f),P H2 (f)} • (4-15) 

This last equation shows that the null channel Jff(/) 
contains less noise power than the difference oinu 1 {f) 1 
nH 2 (f)- The filtering produced by a(f) results in a less 
noisy null estimator than the quantity Uh 1 (/) — nn 2 (/)• 
In the limit that either signal dominates the noise power 
(e.g.,P ffl (/)<P ff2 (/)), 

P z {f) -> max{P ffl (/), P H2 (/)} - min{P ffl (/), P Hs (/)} . 

(4.16) 

In addition, one can form the quantity: 

VpAT) VPhAD + PhAI)' 

As suggested by the label t, this quantity is identical 
to the Student's t statistic, which is used to assess the 
statistical significance of two quantities having different 
means and variances. 



V. CROSS-CORRELATION STATISTICS USING 
COMPOSITE PSEUDO-DETECTOR CHANNELS 

Since the instrumental transcontinental (Ll-Hl, Ll- 
H2) cross-correlations are assumed to be negligible, the 
derivation of the optimal filter when using the pseudo- 
detector channels for Hanford proceeds exactly as has 
been presented in the literature 0, El El with Phi(/)j 
PhAJ) replaced by Pjj(/), P z (f) for the optimal esti- 
mate and the null signal, respectively 



A. Cross-correlation statistic for the optimal 
estimate of the gravitational wave strain 

The cross-correlation statistic is given by 

T/2 



,T/2 

Y Li h ■= / dti I dt 2 s Ll (ti)QL 1 H{ti-t 2 ) S H {t2) : 
J-T/2 J-T/2 

(5.1) 

where T is the observation time and <5i 1 //(t) is the op- 
timal filter, which is chosen to maximize the signal-to- 
noise ratio otY^H- The corresponding frequency domain 
expression is 



Y LlH oc 



Specializing to the case f2 gw (/) = Oo 
mal filter becomes 



Q Ll H(f)=Af Ll H 



dfr L Af)QL lH (f)s H (f). (5.2) 
const, the opti- 

7(1/1) 



\f\ 3 P L Af)PH(f) 



(5.3) 



where MlxH is a (real) overall normalization constant. In 
practice we choose Ml 1 h so that the expected value of the 
cross-correlation is Qq hf 0Q T, where /iioo is the Hubble 
expansion rate Hq in units of £fioo : = 100 km s _1 Mpc -1 . 
For such a choice, 



20n 2 



3H{ 



on 



df 



7 2 (l/l) 



f 6 P L Af)PH(f) 



(5.4) 



Moreover, one can show that the normalization factor 
NliH and theoretical variance, o\ L , of Yl^h are re- 
lated by a simple numerical factor: 



A'' 



LiH 



1 

T 



3^100 

5tt 2 



(5.5) 



1. Limiting case for white coherence and Ph x (/) oc Ph 2 (/) 

If the coherence is white (i.e., PH 1 H 2 {f) = const) and 
the power spectra PhADi Ph 2 (/) are proportional to 
one another, then one can show that the value of the 
cross-correlation statistic Yt i1 h reduces to a linear combi- 
nation of the cross-correlation statistics Yt, 1 h 1 and Yt,^h 2 
calculated separately for Ll-Hl and L1-H2, if we allow 
for instrumental correlations between HI and H2. Thus, 
for this case, combining the point estimates of f2o made 
separately for Ll-Hl and L1-H2 gives the same result as 
performing the coherent pseudo-detector channel analy- 
sis using the single optimal estimator s# (/). 

To show that this is indeed the case, note that 
PHiH 2 {f) = const implies 



T H 1 H 2 (f) ■= \PH l H 2 {f)\" = const. 



(5.6) 
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Representative SI strain sensitivities for 
z = [1 - Ct(f ) - Ot*{f ) + 2|«[f) | 2 ] 1/2 (b (f) - b (f) ) 
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FIG. 4: Same as Fig. Q but for the null signal ZH(f) instead of the optimal estimate Iif(/). Strain spectral densities (i.e., 
absolute value) of 2h(/) (gray or dotted), SHxif) (black), and sjf 2 (/) (dashed), representative of the SI data. Top Panel: 
overlay of the individual amplitude spectral densities with that of the strain spectral density is calculated with the 

SI run-averaged coherence, T s H 1 H 2 {f)- On this scale, the left hand panel shows no discernible difference between the spectra 
for ^H 1 H 2 (f), and with Fh 1 h 2 = 0, suggesting that even the level of coherence seen during the SI run might be sufficiently 
low to allow one to simply combine the Ll-Hl and L2-H2 cross-correlation measurements under the assumption of zero cross- 
correlated noise. The optimality of the estimate 2h(/) is visible here because it is always less than the larger of SH 1 (f) or 
SH 2 (f)- Bottom panel: overlay of individual amplitude spectra with that for |2jj(/)| calculated with ^H 1 H 2 (f) = (i.e., 
assuming no coherence). The difference between the two is very small except for the very lowest frequencies and at narrow line 
features. 



We will drop subscripts for constant quantities. If we 
further assume that Pn 2 {f) = ?7-Ph"i(/)i then 



P Hl H 2 (f) 
PH 2 {f) 



P 



Pk x H 2 (/) 
P Hl (f) 



(5.7) 



Thus, the integrand of the cross-correlation statistic, 

y L Mf)--=slM)QL, H {f)s H {f), (5. 

becomes 
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Y LlH (f) _ iUWlM) [sHAf)(PH 2 (f) - ^ 1 h 2 (/))+^ 2 (/)(Ph 1 (/) - PJW/))] 



LiH 



i -r 



\f?P Ll {f)P Hl {f)PH,{f) (i - r HlH2 (f)) 

Y Ll H 2 (fY 



l-^)%^+(l-A/^) 



where the normalization factor Eq. (|5.4|) is 



3i?. 



1(10 



df 



7 2 (l/l)(^(/) + ffc a (/) - (P Hl H 2 (f) + P* HlH2 (f))) 



pPLAf)P Hl (f)PH 2 (f) (l-T HlHa (J)) 



(i-r) 



Equivalcntly, 



= (i-r) 
= (i-r) 



J y ii«i v r v i> Y Li h 2 



at at 

Y L 1 H 1 IL 1 H 2 



(5.9) 
(5.10) 

(5.11) 
(5.12) 

(5.13) 
(5.14) 



where we used Eq. (|B . f>|> and similar equations to relate NliHi, Nli_h 2 to a 2 LiHi , a 2 LiH2 . 

Substituting the above results for the normalization factors and variances into Eq. (|5.1U|) and integrating over 
frequency, we find: 



Y, 



L X H — 



°Y LlH 

(i-r) 



at at 

*L 1 H 1 * H 2 





(1- 


a\ 







^7 



, P \ Y LlHl „ Y Li h 2 



P 



Y LlHl + a YLiHi (1 - p*^v) Y LlH2 



I 



(5.15) 
(5.16) 

(5.17) 



Or in the notation of Appendix 1X1 



(C 2 2 - Ci 2 ) Yi + (Cn - C 2 i) Y 2 . 

Y L X H - 75 — - -p, ^ , (5.18) 



Cn + C2 



C21 



where Y\ := Yl 1 h 1 , Y2 := Yl 1 h 21 and where we 
used Eqs. (|B1(I . i|B3|) from Appendix [5] to equate 
u y^ Hy i a Y LlH2 with ^lX) C22, and P Hi h 2 /Ph 2 = p/y/rj, 
Ph.hJPh, = P*VV with C12/C22, Cbi/Cu. Thus, we 
see that for the limiting case of white coherence and pro- 
portional power spectra, the pseudo-detector optimal es- 
timator analysis reduces to a relatively simple combina- 
tion of the separate cross-correlation statistic measure- 
ments. 

Finally, note that in the case of zero cross-correlated 



noise (i.e., for p Hl H 2 (f) = 0) we get 

a Y LlH2 Y L lHl + v 2 YliHi Y LlH2 



Yi 



at + at 

a Yl lHl Yl 1 h 1 + a Y ? iH2 Y LlH2 
-2 1 -2 



(5.19) 
(5.20) 



which is the standard method of combining results of 
measurements in the absence of correlations [6J . 



B. Cross-correlation statistic for the null signal 

Once again, the cross-correlation statistic in the fre- 
quency domain is given by 



Y LlZ cx / dfT Li (f)Q LlZ (f)I H (f). 



(5.21) 
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As before, the optimal filter for Sl gw (/) = — const is 1. Limiting case for white coherence and Ph x (/) oc Pn 2 {f) 



Ql i2 (/)=AA LiJ 



7(1/1) 



\f?P L Af)P*U) ' 

where Ml iZ is chosen to be 



20n 2 



df 



7 2 (l/l) 



f 6 P Ll {f)P z {f) 



(5.22) 



(5.23) 



and is related to the theoretical variance oh via: 



T V 5tt 2 



(5.24) 



We start again with the same assumptions that the co- 
herence is white and the power spectra Ph 1 (J), PH 2 (f) 
are proportional to one another (cf. Eqs. Ij5.6|l . (|5. 7\i ). 
Then it is possible to show that the value of the cross- 
correlation statistic Yl iZ reduces to a linear combination 
of the cross-correlation statistics Yt i1 h 1 and Yt, 1 h 2 calcu- 
lated separately for Ll-Hl and L1-H2, if we allow for in- 
strumental correlations between HI and H2. After much 
algebra similar to that presented earlier in Section IV A II 
we obtain: 



Y LlZ 



Jr] (rp - p) 



^L, Ho 



^ 



or, equivalcntly, 



Yl 1z 



\ (v^-p*) (v + v 2 - Vv (p + p*)) V aYL 

Yl 1 h 2 — Yl 1 h 1 



-Vv- 



(5.25) 



(5.26) 
(5.27) 



2. Limiting case for zero cross- correlated noise 



C. Combining triple and double coincident 
measurements 



If also pHiH 2 (f) = 0, then the two interferometer noise 
floors are uncorrclatcd, and the cross-correlation statistic 
Yl iZ for the null channel simplifies further: 



Y, 



Y, 



< v -r <->v 



(5.28) 



Equation (|5.28|l shows that in this limit the quantity 
Yl iZ /&l iZ follows the Student's t distribution. This dis- 
tribution provides a measure to assess the significance of 
the difference between two experimental quantities hav- 
ing different means and variances. Here it provides a 
measure of consistency of the two independent measure- 
ments, Y^^and Yt i1 h 2 : their difference should be con- 
sistent with zero within the combined experimental er- 



In order to make use of this method for the analy- 
sis of future science data, we will need to partition the 
data into three non- overlapping (hence statistically inde- 
pendent) sets: the H1-H2-L1 triple coincident data set, 
and the two Ll-Hl and L1-H2 double coincident data 
sets. The triple coincidence data would be analyzed in 
the manner described in this paper, while the double 
coincidence data (corresponding to measurements from 
different epochs or from different science runs) can be 
simply combined under the assumption of statistical in- 
dependence (cf. Eq. H5.20fl ). 



VI. CONCLUSION 

The approach presented above is fundamentally dif- 
ferent from how the analysis of SI data was conducted 
and represents a manner to maximally exploit the feature 
of LIGO that has two co-located interferometers. This 
technique is possible for the Hanford pair of detectors 
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because, to high accuracy, the gravitational wave signa- 
ture is guaranteed to be identically imprinted on both 
data streams. Coherent, time-domain mixing of the two 
interferometer strain channels can thus be used to opti- 
mal advantage to provide the best possible estimate of 
the gravitational wave strain, and to provide a null chan- 
nel with which any gravitational wave analysis can be 
calibrated for backgrounds. 

An analogous technique of "time-delay intcrfcromctry" 
(TDI) has been proposed in the context of the Laser In- 
terferometer Space Array (LISA) concept 01 El- How- 
ever, in that case the data analysis is very different from 
what is explored in our paper. TDI involves time-shifting 
the 6 data-streams of LISA (2 per arm) appropriately be- 
fore combining them so as to cancel (exactly) the laser- 
frequency noise that dominates other LISA noise sources. 
Even after implementing TDI, the resulting data com- 
binations (with the laser-frequency noise eliminated) are 
not all independent, and may have cross-correlated noises 
from other, non-gravitational-wave, sources. One, there- 
fore, seeks in LISA data analysis an optimal strategy for 
detecting a given signal in these TDI data combinations. 
On the other hand, the method presented in this paper 
is not about canceling specific noise components from 
data; rather, it is about deducing the optimal detection 
strategy in the presence of cross-correlated noise. 

The usefulness of Zuif) is that it may be used to ana- 
lyze cross-correlations for non-gravitational wave signals 
between the Livingston and Hanford sites. This would 
enable a null measurement to be made — i.e., one in which 
gravitational radiation had been effectively "turned off." 
In this sense, using Sif (/) would be analogous to analyz- 
ing the ALLEGRO-LI correlation when the orientation 
of the cryogenic resonant bar detector ALLEGRO is at 
45° with respect to the interferometer arms 0E3. Un- 
der suitable analysis, the cross-correlation statistic Yl iZ 
could be used to establish an "off-source" background 
measurement for the stochastic gravitational wave back- 
ground. 

Ultimately, the usefulness of such a null test will be 
related to how well the relative calibrations between HI 
and H2 are known. If the contributions of h(f) to 'sh 1 (/) 
and sjy 2 (/) are not equal due to calibration uncertain- 
ties, then this error will propagate into the generation 
°f sn(f), 2ff (/). It is possible to estimate this effect 
as follows. Due to the intended use of zh(J) in a null 
measurement, the leakage of h(f) into this channel is the 
greater concern. Considering the structure of Eqs. 1)3. 17[) . 
(|4.3|) . (|4.4|l . it is clear that effects of differential calibra- 
tion errors in sn(f) will tend to average out, whereas 
such errors will be amplified in zn(f). Assume a differ- 
ential calibration error of ±?"(/). Then zn(f) will contain 



a gravitational wave signature 

Sh(f)=2I(f)h(f), (6.1) 

with corresponding power 

5P h {f) = A\Z{f)\ 2 P h {f). (6.2) 
The amplitude leakage affects single-interferometer based 
analyses; the power leakage affects multiple interfer- 
ometer correlations (such as the stochastic background 
search). Assuming reasonably small values for ±7(/), 
if a search sets a threshold on putative gravitational 
wave events detected in channel sjy(/), then the corre- 
sponding contribution in Zn(f) would be approximately 
2|e|p», where |e| denotes the magnitude of the frequency 
integrated differential calibration errors. For any reason- 
able threshold (e.g., p» « 10) above which one would 
claim a detection, and for typical differential calibration 
uncertainties of 2|e| < 20%, then the same event would 
have a signal-to-noise level of p* « 2 in the null channel, 
well below what one would consider meaningful. A more 
careful analysis is needed to quantify these results, since 
calibration uncertainties also propagate into a(/). 

While the focus of this paper is the application of this 
technique to the search for stochastic gravitational waves, 
it appears that any analysis can exploit this approach. It 
should be straightforward to tune pipeline filters and cull 
spurious events by using the null channel to veto events 
seen in the 5i/(/) channel. 
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while the method discussed here combines the data at 
the level of the cross-correlation statistic measurements 
Yt, 1 h 1 an d Yt, 1 h 2 — i.e., after optimal filtering of the in- 
dividual data streams. As such, the method described 
here is not optimal, in general, since it does not take 
advantage of the common gravitational wave strain com- 
ponent h present in HI and H2. However, as shown in the 
main text, when the cross-correlation PH t H 2 {f) is white 
and the power spectra Pjj 1 (/), Pjj 2 (/) are proportional 
to one another, the pseudo-detector optimal estimator 
method reduces to the method described here. 

Consider then a pair of (real-valued) random variables 
Yi, I2 with the same theoretical mean 



A 4 := (Yi) = (Y 2 ) 
and covariancc matrix 



C11 C12 
C21 C22 



where 



{{Y i - l x){Y j - l j)) = {Y i Y j )-^ 



(Al) 



(A2) 



(A3) 



Note that Cyi = C21 since the Yi are real. The absence 
of cross-correlations corresponds to C12 = C21 = 0. 
Now form the weighted average 



Y op t '■— 



(A4) 



It is straightforward to show that Y opt has theoretical 
mean /i opt = (J., and theoretical variance 



(A5) 



Now find the weighting factors Xi that minimize the vari- 



ance of Y opt ■ The result is 



(A6) 



APPENDIX A: GENERAL METHOD OF 
COMBINING MEASUREMENTS ALLOWING 
FOR CROSS-CORRELATIONS 

In this appendix, we present a general method of com- 
bining measurements, allowing for possible correlations 
between them. In the following appendix (Appendix [BJ, 
this method is applied to the case of the Ll-Hl and L1-H2 
cross-correlation statistic measurements, which are taken 
over the same observation period and which may contain 
significant instrumental H1-H2 correlations. 

It is important to emphasize that the method discussed 
in this appendix is not the same as the pseudo-detector 
optimal estimator method discussed in the main text; the 
pseudo-detector method combines the data at the level 
of data streams Sij 1 (/), ^H 2 (f) before optimal filtering, 



or, explicitly, 



Ai = 



C*22 — C\2 

detllCII 



C\\ — C'2 
det||C|| 



(A7) 



where det ||C|| := CnC 2 2 - C 12 C 2 i. 

One can prove the above result by defining an inner 
product 



(A,B) AHCH^B,-, 



* 3 

and rewriting the variance as 



opt 



(C ■ A, C ■ A) 
(C-A,l)2 



(A8) 



(A9) 
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Then t is minimized by choosing A^ such that 
C • A := 2J C l j\ 3 = 1 



(A10) 



for all i. 

For such a choice, 



E. _ Cll + C22 — C*12 — C21 
l ~ fct\\c\\ 

Y opt (C 2 2-C 12 )Y 1 + (C 11 -C 21 )Y 2 



'opt 



det MCI 



so 



y, 



(C 22 — C12) Y\ + (C11 — C21) Y 2 



opt 



C*n + C22 — C12 — C21 
This is the desired combination. 



(All) 
(A12) 

(A13) 



APPENDIX B: APPLICATION OF THE 
GENERAL METHOD TO THE Ll-Hl, L1-H2 
CROSS-CORRELATION STATISTIC 
MEASUREMENTS 

Here we apply the results of the previous appendix 
to the Ll-Hl and L1-H2 cross-correlation measurements. 
We let Y\ denote the cross-correlation statistic Yl 1 h 1 
for the Ll-Hl detector pair, and Y 2 denote the cross- 
correlation statistic Yl 1 h 2 for L1-H2, and assume that 
the measurements are taken over the same observation 
period of duration T. (If the observations were over dif- 
ferent times, then there would be no cross-correlation 
terms and a simple weighted average by a~ 2 would suf- 
fice.) We need only calculate the components of the co- 
variance matrix to apply the method described in the 
previous appendix. 



To calculate the C,j, we assume (as in the main 
text) that the cross-correlated stochastic signal power 
Pn{f) is small compared to the auto-correlated noise in 
the individual detectors, and that there are no broad- 
band transcontinental instrumental or environmental 
correlations — i.e., \P^ 1 h- (/)l ^ s small compared to the 
auto-correlated noise, the cross-correlated stochastic sig- 
nal power, and the H1-H2 cross-correlation \PH 1 H 2 (f)\- 
Then it is fairly straightforward to show that 



L1H1 > 



C22 —^L 1 H 2 



(Bl) 



and 



C 



12 



C- 



21 



C11C22 C11C22 

= }_( 3g ioo N ' 
T V IOtt 2 



df 



l 2 {\.f\)P Hl H 2 {f) 
f 6 P Ll (f)PHAf)PH 2 (f)- 

(B2) 



Note that the above integral is real since Ph 1 h 2 (— /) = 
Ph 1 h 2 (/) an d the integration is over all frequencies (both 
positive and negative). 

If we further consider the limiting case defined by 
white coherence (i.e., pH 1 H 2 {f) = const) and propor- 
tional power spectra (i.e., Ph 1 (/) oc PH 2 (f)), then 
PH lH2 (f)/P Hl (f) and P* HlH2 (f)/PH 2 {f) arc both con- 
stant with values 



P 



11. 



C\ 2 
C 22 



P. 



H!H 2 
P Hl 



C 2 \ 



(B3) 



